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We present an efficient algorithm for generating unitary maps on a d-dimensional Hilbert space 
from a time-dependent Hamiltonian through a combination of stochastic searches and geometric 
construction. The protocol is based on the eigen-decomposition of the map. A unitary matrix 
can be implemented by sequentially mapping each eigenvector to a fiducial state, imprinting the 
eigenphase on that state, and mapping it back to the eigenvector. This requires the design of 
only d state-to-state maps generated by control waveforms that are efficiently found by a gradient 
search with computational resources that scale polynomially in d. In contrast, the complexity of a 
stochastic search for a single waveform that simultaneously acts as desired on all eigenvectors scales 
exponentially in d. We extend this construction to design maps on an n-dimensional subspace of 
the Hilbert space using only n stochastic searches. Additionally, we show how these techniques can 
be used to control atomic spins in the ground electronic hyperfine manifold of alkali atoms in order 
to implement general qudit logic gates as well to perform a simple form of error correction on an 
embedded qubit. 



I. INTRODUCTION 

The goal of quantum control is to implement a non- 
trivial dynamical map on a quantum system as a means 
to achieve a desired task. Historically, the major de- 
velopments in quantum control protocols have been mo- 
tivated by applications in physical chemistry whereby 
shaped laser pulses excite molecular vibrations and rota- 
tions [ , ], and in nuclear magnetic resonance whereby 
shaped rf pulses cause desired spin rotations in magnetic- 
resonance-imaging [ , , ] . More recently, quantum con- 
trol theory has been considered in the development of 
quantum information processors in order to tackle the 
challenges of extreme precision and robustness to noise 
and environmental perturbations [ , , , , Such 
quantum processors are being explored on a wide vari- 
ety of platforms ranging from optics and atomic systems, 
to semiconductors, and superconductors. The design of 
new protocols for quantum control can thus impact a 
wide spectrum applications. 

The simplest approach to quantum control is open- 
loop unitary evolution. In this protocol, the system 
of interest is governed by a Hamiltonian that is a 
functional of a set of time-dependent classical "control 
waveforms", i^[B(t)]. Through an appropriate choice 
of B(t), the goal to reach a desired solution to the 
time-dependent Schrodinger equation at time T, for- 
mally expressed as a time-ordered exponential, U{T) = 

T (^exp |-z H [B(t')] dt'Y^ . The system is said to be 

"controllable" if for any W in the space of unitary maps 
on the Hilbert space of interest, there exists a set of con- 
trol waveforms that such that U (T) = W for some time 
T. Control theorists have long known the conditions on 
H [B(t)] such that the system is controllable in principle, 
but a construction for specifying the desired waveforms is 
generally unknown. The goal of this paper to to provide 



such a construction for a wide class of quantum systems. 
We restrict our attention to Hilbert spaces of finite di- 
mension d. 

Two classes of quantum control problems have been 
primarily considered: state-preparation and full unitary 
maps. In state-preparation, the goal is to map a known 
fiducial initial quantum state to an arbitrary final 
state \i^f). This requires specification of only one col- 
umn of the unitary matrix, i.e. the vector U(T)\ipi)^ as 
compared with the full unitary map, which requires spec- 
ification of all d orthonormal column vectors. The con- 
trast between these tasks is reflected in the complexity 
of numerical searches for the desired waveforms. Opti- 
mal control theory provides a framework for carrying out 
such searches [ ' An objective function J is defined for 
the task at hand, e.g., J[B(t)] = {^l)f\U{T)\^l)i) for state 
preparation or J[B(t)] = Ti: {w^U{T)) for full unitary 
mapping. The optimal controls are the maxima of these 
objective functions. 

In series of papers, Rabitz and coworkers introduced 
the concept of the "control landscape" [11, 12, 13, 14]. 
By discretizing the control functions (e.g., by sampling at 
discrete times), one can treat the objective as a smooth 
function whose domain is a finite set of control variables. 
The topology of this resulting hypersurface governs the 
complexity with which numerical search algorithms can 
find optimal solutions. In the case of state preparation, 
Rabitz et al. showed that for open-loop unitary con- 
trol, the control landscape has an extremely favorable 
topology [ ] . Given a closed-system open- loop Hamilto- 
nian evolution for sufficient time T, all critical points (i.e. 
those values of the control parameters where 5 J = 0) are 
either unit fidelity or zero fidelity; there are neither lo- 
cal optima nor saddle points. Furthermore, the surface 
has a gradual slope as one moves towards the optimal 
points, and there are an infinite number of optimal solu- 
tions connected on a submanifold with large dimension. 
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A^c — 2(i + 2, where Nc > (P — 1 is the number of control 
variables defining the dimension of the overall objective- 
function hyper surface [12]. The lack of false suboptimal 
critical points, the gentle slope, and the flat region near 
a maximum, all enable efficient search algorithms that 
yield fairly robust optimal control waveforms based on a 
simple gradient ascent algorithm from a random seed. A 
collection of random seeds yield a collection of possible 
solutions that can then be further tested for robustness 
to decoherence and noise. 

In contrast, the control landscape for full unitary con- 
trol is less favorable [13]. For SU (d) matrices, there are d 
critical values of the objective function. Of these, there is 
one optimal solution with unit fidelity, an isolated point 
in the control landscape. The remaining d — 1 subopti- 
mal points are saddles. While the lack of local maxima 
may suggest numerical optimization might still be an ef- 
ficient search strategy, empirical studies show otherwise 
[14]. Whereas state preparation search routines converge 
in a number of iterations that is essentially independent 
of (i, the resources necessary for optimization algorithms 
to converge on the full unitary control landscape grow ex- 
ponentially with d. Brute force search is thus a very poor 
strategy for full unitary control on all but the smallest 
dimensional Hilbert spaces. 

Explicit constructions for full unitary control have 
been established in special cases where the form of the 
Hamiltonian allows. For example, Khaneja et al. showed 
that the problem of generating unitary matrices on a sys- 
tem of weakly coupled qubits can be reduced to the so- 
lution of a geodesic equation [ ]. Brennen et al. showed 
that by considering controls on overlapping 2-d subspaces 
it is possible to create arbitrary controls through Givens 
rotations [ ] . Such constructive procedures are less com- 
putationally intensive than their random search counter- 
parts, and moreover, yield control fields that are more 
physically intuitive. They are, however, restricted to 
control systems with particular structures and are not 
applicable in more generic cases. 

In this paper we develop a hybrid protocol for full uni- 
tary control that combines efficient numerical search pro- 
cedures with constructive algorithms, applicable for any 
finite dimensional Hilbert space with minor restrictions, 
thereby extending the work of Luy et al. [ ] . We lever- 
age off of the efficiency of numerical searches for wave- 
forms that generate a desired state mapping. Our proce- 
dure requires only d stochastic searches and the length of 
the resulting control sequence is approximately 2d times 
the time of a single state preparation. Our work is anal- 
ogous to that of [ ], where state-preparation provides 
a good starting point for iterative searches. Our pro- 
cedure, however, is fully constructive and deterministic 
once appropriate state mappings are found. 

The remainder of this paper is structured as follows. 
In Sec. II we present our hybrid protocol for construct- 
ing general unitary maps by combining efficient numeri- 
cal searches with a deterministic algorithm. In addition 
to unitary maps on the full Hilbert space, this scheme al- 



lows us to construct maps on a subspace with a complex- 
ity that scales as the dimension of that space. Finally, 
in Sec. HI, we apply our unitary matrix construction to 
control the large manifold of magnetic sublevels in the 
ground electric states of an alkali atom (e.g. ^^^Cs) [17]. 
We show how to construct a set of unitary matrices on 
SU{d) that are often considered as qudit logic gates in 
a fault-tolerant protocol. In addition, we apply our con- 
struction for subspace mapping to encode logical qubits 
in our qudit, and simulate an error correcting code that 
protects against magnetic field fluctuations. 

II. UNITARY CONSTRUCTION 

In this section we define an efficient protocol for con- 
structing arbitrary unitary maps based on state prepara- 
tion. Any unitary matrix has an eigen-decomposition, 

U = J2e-^^^\^j){cPj\ = J]e-^^^l*^-><*^l, (1) 

j 3 

where in the second form we expressed /7 as a product 
of commuting unitary evolutions by moving the projec- 
tors into the exponential. A general unitary map can 
be thus be constructed from d propagators of the form 
exp{— zAj|0j)((^j|}, one for each eigenvalue/eigenvector 
pair. These unitary propagators can now be constructed 
using state mappings. We begin by noting that there 
exists some Vj G SU (d) that satisfies 

g-iA,-|0,-)(0,-| ^ ^-iXjVj\0){0\Vj ^ l/tg-iA,-|0)(0|y^.^ 

where |0) is a fixed "fiducial state". The sole require- 
ment on Vj is that KO|V^|0j)p = 1, i.e., it must be a 
mapping from to |0). Therefore, we can create the 
unitary propagator exp{— zAj (^j |} by using a state 
preparation to map the eigenvector of [/, onto the 
fiducial state |0), applying the correct phase shift, and 
finally mapping the fiducial state back to the eigenvec- 
tor with the time-reversed state preparation. A general 
unitary map is thus constructed via the sequence, 

U = Vje-'^'\'^^'\Vd . . . F2^e-^^^|o>^oiF2^i^e-^^^|o>^oiFi. 

(3) 

Each of the propagators Vj is specified by a control wave- 
form that generates a desired state mapping. One can 
efficiently find such control fields based on a numerical 
search that employs a simple gradient search algorithm, 
as described above. To generate an arbitrary element of 
SU{d)^ we require at most d such searches. Moreover, 
the full construction consists of 2d state preparations in- 
terleaved with d applications of the phase Hamiltonian, 
leading to an evolution that is only of order d times longer 
than a state mapping evolution. 

This construction places only two requirements on the 
Hamiltonian in addition to controllability. Firstly, the 
dynamics must be reversible such that if we can gener- 
ate the unitary evolution V^-, we can trivially generate 
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the unitary by time-reversing the control fields. Note 
that this is not the same as finding a state preparation 
that goes in the opposite direction, |0) there are 

many unitary propagators that map |0) so it is 

unlikely to find the unique operator Vj from a stochastic 
search. Secondly, we require access to a control Hamil- 
tonian that applies an arbitrary phase to one particular 
fiducial state |0) relative to all of the remaining states in 
the Hilbert space, exp{— zAj |0)(0|}. This latter require- 
ment is the most restrictive, but can be implemented in 
a wide variety of systems. An example is discussed in 
Sec. III. 



A. Subspace Maps 

We have so far considered two kinds of maps on our 
d-dimensional Hilbert space H: d x d unitary matrices 
and state-to-state maps. The former corresponds to a 
map U : ^ 7Y, while the latter specifies a map on a 
one-dimensional space. Intermediate cases are also im- 
portant. In particular, we are often interested in unitary 
maps that take subspace A of arbitrary dimension n to 
subspace S, according to T : A ^ B. Examples include 
the encoding of a logical qubit into a large dimensional 
Hilbert space {A 7^ B) and a logical gate on encoded 
quantum information {A = B). Above we showed that 
the design of a fully specified unitary matrix required 
search for d waveforms that define d state preparations 
(trivially a state mapping requires one such search). We 
show here how unitary maps on subspaces of dimension 
n can be constructed from exactly n such numerical so- 
lutions. 

Formally, a unitary map between two subspaces A and 
B of dimension n is defined as a map between between 
their orthonormal bases {\ai)} and 

n 

T^{A^B) = J2\bi){ai\(BU±, (4) 

where U± is an arbitrary unitary map on the orthogonal 
complement A± whose dimension is d — n. State prepa- 
ration is the case n = 1; a full unitary matrix is specified 
when n = d. Clearly for n d the map is not unique, 
with implications for the control landscape and the sim- 
plicity of numerical searches described above. As a first 
naive construction of T{A S), one might consider a 
sequence of one-dimensional state mappings, 

n 

T^{A^B)ll[Ti{\ai)^\bi)). (5) 

This does not, however, yield the desired subspace map 
because each state mapping acts also on the orthogonal 
complement, so, e.g. \bi) is affected by Ti (|a2) ^ I62)) 
and subsequent maps will move formerly correct basis 
vectors to arbitrary vectors in the orthogonal compo- 
nent. We can resolve this problem by instead construct- 



ing subspace maps as a series all well-chosen rotations 
that maintain proper orthogonality conditions. 

To construct the necessary unitary operators, we make 
use of the tools described above: arbitrary state mapping 
based on an efficient waveform optimization and phase 
imprinting on a fiducial state. With these, we define the 
unitary map between unit vectors \a) and 

5(|a),|6))^e--l*H0l=j_2|0)(0|. (6) 

Here = 7V(|a) — where we have chosen the 
phases such that {b\a) is real and positive, and = 
2 (1 — {b\a)) is the normalization. This unitary operator 
has the following interpretation. In the two-dimensional 
subspace spanned by \a) and S is a 7r-rotation that 
maps S\a) = \b). In contrast to the state preparation 
map, Eq. (4) with n = 1, this map acts as the identity on 
the orthogonal complement to the space. This property 
is critical for the desired application. 

With these 2D primitives in hand, we can construct 
the subspace map according to the prescription, 

Tn{A^B)=Sn-.-S2Si, (7) 

where Sk = S {\dk), \bk)) and 

\dj) =Sj-i...S2Si\aj). (8) 

This sequence does the job because each successive rota- 
tion leaves previously mapped basis vectors unchanged. 
To see this, we must show that at step j, the basis vectors 
{1^1)7 1^2)7 • • • 7 cire unchanged by Sj. This will be 

true when this set is orthogonal to the vectors \dj) and 
\bj). Orthogonality to | bj ) is trivial since the basis vectors 
of B are orthonormal. We must thus prove, {dj\bk) = 0, 
Vj > k. We can do this by induction. For an arbi- 
trary /c, assume the conjecture is true for all j such that 
jo > j > k, and thus Sj\bk) = \bk) up to j = jo. This 
implies that {djQ-^i\bk) = since, 

= («jo+ii4---4i^fc) 

= (%o+il«/c) =0. (9) 

To complete our proof by induction, we must show that 
for any /c, the conjecture is true for j = k -\- 1. This 
follows since, 

(a/e+i|6/c) = {ak ^ l\s\sl . . . sl\bk) 

= {ak+i\ak)=0. (10) 

With this protocol we can construct unitary maps on a 
subspace of dimension n with optimized waveforms that 
corresponded to exactly n prescribed state preparations. 
In the following section we apply these tools to qudit 
manipulations in atomic systems. 
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FIG. 1: The hyperfine structure of ^^^Cs in the 6S1/2 ground 
state. Microwaves (blue) and rf magnetic fields (purple) 
provide controllable dynamics on the 16-dimensional Hilbert 
space. A detuned laser light shift (red) can be used to create 
a relative phase between the F — A and F — manifolds. By 
considering controls on the subspace of the orange states we 
recover a system that satisfies the criteria proposed in Sec. II. 



In addition to an efficient method for designing and im- 
plementing state-to-state mappings, our protocol places 
certain requirements on the available control tools. 
Firstly, the system dynamics must be reversible so that 
we can trivially invert a state mapping. This is easily 
achieved through phase control. Secondly, we require 
phase imprinting on a single fiducial state. While this 
cannot be accomplished using solely microwave and rf- 
control, by introducing an excited electronic manifold, 
an off-resonant laser-induced light-shift can achieve this 
goal. We restrict our system to one spin manifold (here 
the F = 3, but in principle either will do) and a single 
state from F = 4 manifold, e.g. |F = 4, m = 4), which 
acts as the fiducial state. By detuning far compared to 
the excited state line width of 5 MHz, but close com- 
pared to the ground- state hyperfine splitting of 10 GHz, 
we imprint a light shift solely on the |F = 4, m = 4) state 
with negligible decoherence. Using rf- magnetic fields to 
perform rotations in the F = 3 manifold, and microwaves 
to couple to the fiducial state, we obtain controllable and 
reversible dynamics. Note that we may include the fidu- 
cial state in our Hilbert space, for a total of 8 sublevels, 
or treat it solely as an auxiliary state and restrict the 
Hilbert space to the 7-dimensional F = 3 manifold. 



III. APPLICATIONS TO ATOMIC SPIN 
CONTROL 

In this section, we apply our results to the control 
of the ground-electronic manifold of magnetic sublevels 
in alkali atoms. Atomic spins are natural carriers of 
quantum coherence for use in various quantum informa- 
tion processing applications [18, 19, 20, 21]. In pre- 
vious work we showed that the full ground-electronic 
subspace of coupled electron and nuclear spin can be 
rapidly controlled through combinations of static, radio, 
and microwave-frequency ac-magentic fields, with negli- 
gible decoherence [ ]. A schematic for the specific case 
of ^^^Cs, with nuclear spin / = 7/2 and two ground- 
electronic hyperfine manifolds with total angular momen- 
tum F = 3 and F = 4, is shown in Fig. 1. A static 
bias magnetic field breaks the degeneracy and specifies 
an rf-resonance frequency by the Zeeman splitting in a 
given manifold. Control of the amplitude and phase of 
the rf-magnetic fields oscillating in two spatial directions 
allows one to independently rotate these two manifolds. 
Resonant microwaves can be used to excite transitions 
between F = 3 and F = 4, driving coherent SU{2) ro- 
tations between two magnetic sublevels, as specified by 
a given (nondegenerate) transition frequency. Such con- 
trols together can be used to generate an arbitrary uni- 
tary transformation on the d — 2(2/ + 1) = 16 dimen- 
sional Hilbert space. In our previous work we showed how 
we could design state-preparation mappings through sim- 
ple gradient searches [17]. In the present work we show 
how we can employ this tool to design more general uni- 
tary maps based on the protocol of Sec. II. 



A. Constructing qudit unitary gates 

The standard paradigm for quantum information em- 
ploys two-level systems - qubits - in order to implement 
binary quantum-logic based on SU{2) transformations. 
Extensions beyond binary encodings in d > 2 system 
- qudits - based of SU{d) transformations have also 
been studied and may yield advantages in some circum- 
stances [io, z^, 2o]. Of particular importance for fault- 
tolerant operation is implementation of these transforma- 
tions through a finite set of "universal gates". Our goal 
here is to show how important members of the universal 
gate set can be implemented using our protocol. 

In choosing a universal gate set appropriate for error 
correction, it is natural to consider generalizations of the 
Pauli matrices X and Z which generate SU (2). The gen- 
eralized discrete Pauli operators for SU (d) are defined 

x\j) = U©i) 

Z\j) = (11) 

Here refers to addition modulo d and uo is the primitive 
dth root of unity, uj = ex.p{i27r / d} . By considering the 
commutation relation of X and Z, the remaining gen- 
eralized Pauli operators have the form u^X^ ^ defining 
the elements of Pauli group for one qudit (up to a phase) . 
This group is a discrete (finite dimensional) generaliza- 
tion of the Weyl-Heisenberg group of displacements on 
phase space. 

Another important group of unitary matrices in the 
theory of quantum error correction is the Clifford group, 
given its relationship to stabilizer codes [ ]. These group 
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elements map the Pauli group back to itself under conju- 
gation. Expressed in terms of their conjugacy action on 
X and Z, the generators of the Clifford group for single 
qudits are 



HXH^ = Z, 
SXS"^ = xz, 
GaXGi=X\ 



HZH^ = X-^ (12) 

SZS"^ = Z (13) 
GaZGi = Z-" 

when gcd(a,(i) = 1 (14) 



H and S are direct generalization of the Haddamard and 
phase-gates familiar for qubits [ ]. The (i-dimensional 
H is the discrete Fourier transform 



Vd ^ 

and 5 is a nonlinear phase gate 

S\j) = J- odd, 

S\j) = uj^ /"^IJ) j even. 



(15) 



(16) 
(17) 



The operator Ga is a scalar multiplication operator with 
no analog in the standard Clifford group on qubits, de- 
fined by 



Ga\j) = 



(18) 



where the multiplication is modulo d. The only such 
multiplication operator for 2-level systems is the identity 
operator. 

While both the generalized Pauli and Clifford groups 
have utility in quantum computing, it is clear from their 
descriptions that unlike their qubit SU (2) counterparts, 
these unitary matrices do not arise naturally as the time 
evolution operators governed by typical Hamiltonians. 
This fact is not relevant to our unitary construction, 
which requires only knowledge of the operators' eigenvec- 
tors and eigenvalues. Using the time-dependent Hamilto- 
nian dynamics with couplings illustrated in Fig.l we have 
engineered control fields to create the generators of both 
the Pauli and Clifford groups acting on the 7-dimension 
F = 3 hyperfine manifold. The duration of waveforms 
is approximately 1.5 ms, which is significantly shorter 
than the decoherence time of the system. In principle, 
the durations of these waveforms could be decreased by 
an order of magnitude or more by using more powerful 
control fields. Our objective function for creating a de- 
sired unitary W is the trace distance J[W] = Tr (W^U)^ 
where U is the unitary matrix generated by our con- 
trol waveforms. Based on our protocol, employing state 
mappings that have fidelities of 0.99, our construction 
yields unitary maps that reach their targets with fideli- 
ties of J[Z] = 0.9866, J[X] = 0.9872, J[H] = 0.9854, 
J[S] = 0.9892 and JlGs] = 0.9801. 

As an example, in Fig. 2 we show the control sequence 
for the discrete Fourier transform. The unitary map gen- 
erated by this sequence should act to transform eigen- 
states of Z into eigenstates of X and vice versa. We 
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FIG. 2: Optimized control fields for implementing the 7- 
dimensional Fourier transform on the F — 3 hyperfine man- 
ifold in ^^^Cs. The duration of the pulse is less than 1.2 ms 
and yields a unitary map that has an overlap of 0.9854 with 
the desired target. As an example, we show the action of the 
resulting unitary on the Z-eigenstates of angular momentum. 
The conjugate variable of Fz is the azimuthal angle (f). If we 
Fourier transform a Z-eigenstate, a state with a well defined 
value of Fz , we obtain a state that has a well defined value of 
0, a squeezed state. 



illustrate this through a Wigner function representation 
on sphere [ ]. The Z eigenstates are the standard basis 
of magnetic sublevels, whose Wigner functions are con- 
centrated at discrete latitudes on the sphere. Fig. 2a. Ap- 
plying the control fields to each of these states yields the 
conjugate states, with Wigner functions shown in Fig. 2b. 
These have the expected form. They are spin squeezed 
states concentrated at discrete longitudes conjugate to 
the Z eigenstates. The Z and X eigenstates are analo- 
gous to the number and phase eigenstates of the harmonic 
oscillator in infinite dimensions. 
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B. Error-correcting a qubit embedded in a qudit 

The ability to generate unitary transformations on 
two-dimensional subspaces allows us to encode and ma- 
nipulate a qubit in a higher dimensional Hilbert space 
in order to protect it from errors. Such protection can 
take a passive form through the choice of a decoherence- 
free subspace [26, 27], or active error correction through 
an encoding in a logical subspace chosen to allow for 
syndrome diagnosis and reversal [28, 29]. Typically, er- 
ror protection schemes involve multiple subsystems (e.g. 
multiple physical qubits) to provide the logical subspace. 
While tensor product Hilbert spaces are generally neces- 
sary to correct for all errors under reasonable noise mod- 
els, for a limited error model, one can protect a qubit by 
encoding it an a higher dimensional qudit [ ]. We con- 
sider such a protocol as an illustration of our subspace- 
mapping procedure. 

As an example, we consider encoding a qubit in the 
ground-electronic hyperfine manifold of ^^^Cs and pro- 
tecting it from dephasing due to fluctuations in external 
magnetic fields. In the presence of a strong bias in the 
z-direction, the spins are most sensitive to fluctuations 
along that axis. For hyperflne qubits, one solution is to 
choose the bias such that two magnetic sublevels see no 
Zeeman shift to first order in the field strength (a "clock 
transition"). An alternative is to employ an active error 
correction protocol analogous to the familiar phase-fiip 
code [ ]. 

We take our "physical qubit" computational basis to 
be the stretched states, |0) = |3,32) and |1) = 14,4^), 
states easily prepared via optical pumping and controlled 
via microwave- drive rotations on the Bloch sphere. Here 
we have used the shorthand labeling the two quantum 
numbers |F, m^), and have denoted the relevant quan- 
tization axis by the subscript on the magnetic sublevel. 
Such states, however, are very sensitive to dephasing by 
fiuctuations along the bias magnetic field, and such er- 
rors are not correctable. We choose as our encoded qubit 
basis stretched states along a quantization axis perpen- 
dicular to the bias (x-axis), {|0) = |3, 3a^), |1) = |3, — 3^;)}. 
Choosing this basis, a dephasing error in the z-direction 
acts to transfer probability amplitude into an orthogonal 
subspace. Such errors that can be detected and reversed 
without loss of coherence. 

Our error correction protocol works as follows (see 
Fig. 3). Consider an encoded qubit = a\0) + /3\1). 
The error operator due to B-field fiuctuations is the gen- 
erator of rotations, F^. Assuming a small rotation angle 
2e, when such an error occurs, the state of our encoded 
qubit is mapped to 

IVi) « 1^/5) + e (a|3, 2,) + /3|3, -2,)) . (19) 

The error acts to spread our qubit between two orthog- 
onal subspaces, \mx\ = 3 and \mx\ = 2. To diagnose the 
syndrome we must measure the subspace without mea- 
suring qubit. We can achieve this by coherently map- 
ping the error subspace to the upper hyperfine mani- 



fold, followed by a measurement that distinguishes the 
two hyperfine manifolds, F = 3 and F = 4. Such 
a coherent mapping cannot be achieved through sim- 
ple microwave- driven transitions since the bias field is 
along the z-direction while the encoded states are mag- 
netic sublevels along the x-direction. We can instead 
use the construction of unitary operators on a subspace 
described in Sec. H to design 7r-rotations that take the 
error states to the upper manifold. This is tricky for our 
implementation because our protocol only included one 
magnetic sublevel in the F = 4 manifold so as to en- 
sure proper phase imprinting. The solution is to switch 
the auxiliary state in the upper manifold between two 
different subspace maps. First, we consider the control 
system where 14,4;^) is our auxiliary state and perform a 
TT-rotation that maps \3^2x) to 14,4^), leaving the rest of 
the space invariant. Then employ control on the system 
where |4, — 4^) is the auxiliary state and map |3, —2x) 
to 1 4, — 4^), with the identity on the remaining space. A 
QND measurement of F collapses the state to the initially 
encoded state when the measurement result is F = 3, or 
to the state a|4,4a,) + /3|4, -4^,), if we find F = 4. In the 
final step of the protocol, if an error occurred, we con- 
ditionally move the error subspace back to the encoded 
subspace, which can be achieved through reverse maps 
of the sort described above. 

We simulate here the coherent steps in the error cor- 
rection protocol. These are implemented through our 
efficient search technique to construct subspace maps for 
the sequences 

{|4,4,),|3,3,)} ^ {|3,3,),|3,-3,)} 
{|3,2,),|3,-2,)} ^ {|4,4,),|4,-4,)} 
{|4,4,),|4,-4,)} ^ {|3,3,)|3,-3,)} 

Each of these maps are achieved through a sequence of 
SU{2) TT-rotations on a two-dimensional subspace that 
leave the orthogonal subspaces invariant. Starting from 
numerical searches for state preparation maps that have 
fidelity greater than 0.99, we obtain subspace maps with 
comparable fidelities. The performance of this error cor- 
rection procedure is shown in Fig. 3B. We plot the fidelity 
between the initial state and the post-error-corrected 
state, averaged over random initial pure states of the 
physical qubit, versus the magnitude of the error as de- 
scribed by the rotation angle induced the stray magnetic 
field. Even with imperfect subspace transformations the 
error correction protocol is significantly more robust than 
free evolution. Of course, like all quantum error correc- 
tion protocols, we assume here that the time necessary 
for diagnosing the syndrome and correcting an error is 
sufficiently shorter than the dephasing time, so that the 
implementation of error correction does not increase the 
error probability. 

In practice, the most challenging step in the error cor- 
rection protocol in this atomic physics example is mea- 
surement of the syndrome. This requires addressing of 
individual atoms and measuring the F quantum num- 
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FIG. 3: (A) A schematic of the error correction proto- 
col we have designed using subspace maps. We track 
the basis elements of our encoded subspace, here |0) is 
red and |1) is blue, via their populations in the x and z 
bases, (i) The initial embedded qubit we wish to protect 
is in a superposition of the 14,4^) and \3,3z) states, (ii) 
We use subspace maps to encode the state in the basis 
|3, 3a;),|3, — 3cc). (iii) In this basis a small z-rotation shifts 
population into the states |3,2x) and |3, — 2^). (iv) Us- 
ing subspace maps we can transfer the small population 
that has left the encoded space in the to states, 14,4^) 
and |4, — 4^). Now we can perform a non-demolition mea- 
surement of the total angular momentum F. If F = 3 
we can be certain our state lies in the encoded subspace 
(ii). If we measure F = 4, the system is in configuration 
(v), which we then conditionally transform back to the 
encoded state. In (B) we examine the performance of the 
error correction. On the x-axis we have the angle of ro- 
tation in the ^-direction due to the magnetic field error. 
On the ?/-axis is the fidelity between the initial and post 
error states, as average over pure states drawn from the 
Harr measure. The blue line shows the fidelity of the er- 
ror corrected states and the green the fidelity if the state 
had simply stayed in the subspace 14,4^), 13,3^). 



ber in a manner that preserves coherence between mag- 
netic sublevels. In principle, this can be achieved through 
a QND dispersive coupling between an atom and cav- 
ity mode that induces an F-dependent phase shift on 
the field that could be detected [ ]. Alternatively, F- 
dependent fluorescence from a given atom would allow 
this code to perform "error detection", without correc- 
tion. 



IV. SUMMARY AND OUTLOOK 



In this paper we have presented a protocol for con- 
structing unitary operators that combines the strengths 
of both stochastic and geometric control techniques. By 
utilizing stochastic searches to construct state prepara- 
tions, as opposed to stochastically searching for full uni- 
tary maps, our protocol requires computational resources 
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that scale only polynomial with the dimension of the 
Hilbert space of our system. The length of the control 
pulses also scales polynomially with d. These stochastic 
search techniques place only very mild restrictions on the 
types of Hamiltonian controls with which our protocol is 
applicable. Additionally, the controls easily generalize to 
the case where one wishes to control only a subspace of 
a larger Hilbert space. For subspace control, the num- 
ber of searches required scales as the dimension of the 
subspace, not as that of the embedding Hilbert space. 

Hybrid stochastic/geometric control schemes yield a 
very promising path towards unitary control sequences 
that balance broad applicability with ease of implemen- 
tation [ ]. The most restrictive element of our protocol 
is the requirement that we can impart a desired phase 
on one and only one state (a U{2) operation). A much 
less restrictive procedure is to employ a control Hamil- 
tonian that acts in imprint a relative phase between two 
states in the Hilbert space (an SU(2) operation). This 
type of operation could be implemented through, e.g., 
the microwave controls described in Sec. HI. As a gen- 
eralization of our protocol, eigenstates would be mapped 
pairwise to two chosen fiducial states where an external 
field generates the desired phase difference. The difficulty 
with this approach is that we require stochastic searches 



for a control waveforms that maps a 2D subspace in 
one step, rather than than our two-step procedure which 
maps each basis vector separately. The topology of the 
control landscape for such waveforms and complexity of 
such a stochastic search are not known, though we expect 
this to be polynomial in d. 

While we have primarily emphasized here an exponen- 
tial speedup in the search for control waveforms that gen- 
erate unitary maps, a constructive protocol brings addi- 
tional possible advantages. By exploiting the geometry of 
a problem, we can engineer robustness more easily than 
in a stochastic setting. For example, the microwave and 
rf-controls discussed in Sec. HI consist of representations 
of SU{2) rotations in different subspaces of the Hilbert 
space. There are well known composite pulse techniques 
that implement rotations SU (2) that are robust to errors 
in the individual pulse amplitudes and detunings [ , 32] . 
In future work we will explore protocols that import these 
methods in order to efficiently search for and implement 
robust SU{d) transformations. 

This research was supported by NSF Grants No. 
PHY-0653599 and No. PHY-0653631, ONR Grant No. 
N00014-05-1-420, and lARPA Grant No. DAAD19-13- 
R-0011. 
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